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Abstract. We solve the Riemann problem for the deceleration of arbitrarily magnetized relativistic ejecta injected into a static 
unmagnetized medium. We find that for the same initial Lorentz factor, the reverse shock becomes progressively weaker with 
increasing magnetization ct (the Poynting-to-kinetic energy flux ratio), and the shock becomes a rarefaction wave when <J 
exceeds a critical value, <J C , defined by the balance between the magnetic pressure in the ejecta and the thermal pressure in 
the forward shock. In the rarefaction wave regime, we find that the rarefied region is accelerated to a Lorentz factor that is 
significantly larger than the initial value. This acceleration mechanism is due to the strong magnetic pressure in the ejecta. 
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INTRODUCTION 

Relativistic jets are believed to exist in active galactic nuclei (AGNs), black hole binaries, and gamma-ray bursts 
(GRBs), but their composition is still poorly understood. It has been argued that magnetic fields could play an important 
dynamic role in these jets [1, 2, 3, 4], but the degree of magnetization, quantified by the magnetization parameter cr 
(the ratio of electromagnetic to kinetic energy flux), is poorly constrained by observations. GRB afterglow modeling 
indicates that the ejecta are more magnetized than the ambient medium, suggesting a possibly important dynamic role 
for magnetic fields in GRB jets [5, 6, 7, 8] . A useful diagnostic for the degree of jet-magnetization can be obtained from 
the interaction between the decelerating jet and the ambient medium. Added magnetic field pressure in the jet alters 
the condition for formation of a reverse shock (RS) as well as its strength [9]. Analytical studies of the deceleration 
of a GRB fireball with arbitrary magnetization [10] suggest novel behavior that does not exist in pure hydrodynamic 
(HD) (a = 0) models [11, 12], However, consensus on the conditions required for the existence of the RS or how 
Poynting flux is transferred to kinetic flux in the interaction region has not yet been achieved [10, 13, 14]. We present 
a one-dimensional study of the interaction between a magnetized relativistic flow and a static, unmagnetized external 
medium. A Riemann problem is solved analytically over a broad range of cr. 


RIEMANN PROBLEM 

We consider a Riemann problem consisting of two uniform initial states (left and right) with discontinuous hydro- 
dynamic properties specified by the rest-mass density p, gas pressure p, specific internal energy u, specific enthalpy 
h = 1 + u/ pc 2 + p/ pc 1 , and normal velocity v N . The right state (the medium external to the ejecta) is assumed to be 
a cold fluid with constant density, at rest. Specifically, we select the initial conditions: p R = 1.0p 0 , p R = 1 0 2 p 0 c 2 , 
v r = v r = 0.0, where p 0 is an arbitrary normalization constant (our simulations are scale-free) and c is the speed of 
light. The left state (the propagating relativistic ejecta) is assumed to have a higher density and pressure than the right 
state, as well as a relativistic velocity. Specifically, p L = 10 2 p 0 , p L = 1.0p 0 c 2 , and = v£ = 0.995c ( y L ~ 10). The 
fluid is described by an adiabatic equation of state p p r with T = 4/3. 

To investigate the effects of magnetic fields, we consider a perpendicular field component in the ejecta with 
B y = 31.623, 100.0, 316.23, and 447.21 in units of (47TPQC 2 ) 1 / 2 measured in the laboratory frame, corresponding to 
<7 = B 1 /Any 1 hpc 1 ~ B 1 /Any 1 pc 1 being 0.1, 1.0, 10.0, and 20.0, respectively. This field is motivated by the predicted 



toroidal field domination at the deceleration radius for GRB outflows [10]. Increasing cr increases the total (kinetic 
plus magnetic) energy density of the left (ejecta) state. 


RESULTS 

We calculate exact solutions of this problem, using the code of Giacomazzo & Rezzolla (2006) [15], in the region 
0.8 < x < 1.2 with an initial discontinuity at x = 1 .0, where x is in arbitrary units [16]. 

For a = 0.1, the solution shows a right-moving fast shock (FS: forward shock; S _>), a left-moving fast shock (RS: 
reverse shock; ^S) relative to the contact discontinuity (C). For cr = 1.0, the solution shows similar profiles (, SCS _>) 
as for cr = 0.1. The FS is stronger (due to a higher jump in pressure) and slower (more deceleration relative to the 
frame of the contact discontinuity), while the RS is weaker but faster. These features are expected from analytical 
work [10, 14], and agree with 1-D relativistic MHD simulations [17, 18]. When the magnetization of the flow exceeds 
cr = 2.7, the shock profiles change drastically. For cr = 10.0 and cr = 20.0, a prominent left-going rarefaction wave 
L R) is observed, instead of a left-going shock. When the rarefaction wave propagates into the ejecta, density and 
gas pressure decrease, and the flow velocity increases. The terminal Lorentz factor of the left (ejecta) state and the 
FS region reaches y ~ 14 for cr = 10 and y > 16 for cr = 20. This magnetic acceleration mechanism stems from the 
magnetic pressure in the ejecta 1 . 

This magnetic acceleration mechanism is solely an MHD effect and a strong magnetic field is required to generate 
the rarefaction wave. This is different from the HD/MHD boost mechanism [20, 21]. 

In the reverse shock cases (cr = 0.1, 1), the upstream magnetic pressure is lower than the gas pressure in the forward 
shock, while in the rarefaction wave cases (cr = 10, 20), the upstream magnetic pressure exceeds the gas pressure in the 
FS. Thus, the balance between the upstream magnetic pressure in the unshocked ejecta region and the FS gas pressure 
in the shocked medium [10, 19] provides the condition separating the two regimes. This condition can be derived 
analytically. For the interaction between relativistic ejecta and an external medium, there exist four physically distinct 
regions: (1) unshocked medium, (2) shocked medium, (3) shocked ejecta, and (4) unshocked ejecta. A critical a c value 
is given as O c = 2p l (y 4 — 1) (4y 4 + 3)/3p 4 — 8 Pj y| /3 p 4 . The condition for the existence of a reverse shock is <7 < a c . 
The condition for a rarefaction wave and magnetic acceleration is a > a c . We adopted p l = p R = 1 .0, p 4 - p L = 10 2 , 
and y 4 = y L = 10.0, so that the critical value is a c ~ 2.7. Our calculations indicate that a c marks the transition point 
where neither a reverse shock nor a rarefaction wave is established. 

To better understand the magnetic acceleration mechanism, we plot the Lorentz factor as a function of a in Fig. 1(a). 
For the magnetic acceleration case, this is the terminal Lorentz factor after acceleration. Because of the dependence 
of o c on y L , a higher cr is needed to achieve acceleration for a higher y L . The terminal Lorentz factor can be estimated 
analytically by requiring that the thermal pressure in the FS region balance the magnetic pressure in the region 
through which the rarefaction wave has propagated. For the terminal Lorentz factor y , this condition can be expressed 
roughly as y ~ (3y|crp 4 /8p | ) 1 4 . Crosses in Fig.l(a) denote values of estimated terminal Lorentz factors for model 
parameters, y 4 = y L = 20, p , = p R = 1.0, and p 4 = p L = 10 2 , in good agreement with the exact solution of the Riemann 
problem in the reverse rarefaction wave regime. 

To investigate the acceleration efficiency, we present in Fig. 1(b) the terminal Lorentz factor y, and its ratio to the 
initial Lorentz factor (y/y L ) as a function of the initial ejecta Lorentz factor y L . While a ejecta with a higher initial 
Lorentz factor reaches a higher terminal Lorentz factor, a lower initial Lorentz factor implies a higher acceleration 
efficiency. From the equation for the terminal Lorentz factor, it follows that y/y 4 — (3cp 4 /8p 1 ) 1/,4 y ) j 1 ' /2 , in good 
agreement with the exact solution of the Riemann problem in relativistic regime. 


DISCUSSION 

Our results have implications for understanding deceleration of strongly magnetized outflows, possibly present in 
GRBs and AGNs. Exact solutions indicate that the condition for the existence of a reverse shock is a < a c [10, 14]. 
The paucity of bright optical flashes in GRBs [22] may, among other interpretions, be attributed to highly magnetized 
flows. Furthermore, the magnetic acceleration mechanism discussed here suggests that a and y are not independent 
parameters at the deceleration radius. For high - a flows, the ejecta would experience magnetic acceleration at small 


1 We note that Romero et al. (2005) [19] also discovered the rarefaction wave regime discussed in this paper, but did not investigate the magnetic 

acceleration mechanism and its astrophysical implications in detail. 





FIGURE 1 . (a) The c-dependence of the maximum Lorentz factor in the shocked region in the case with initial Lorentz factor 

y L = 20. (b) The dependence of the acceleration efficiency y t /y L on the initial Lorentz factor y L with initial magnetization o = 100. 
Solid lines refer to the RR regime and dotted lines to the RS regime obtained with the exact solution. Crosses are the values of the 
estimated terminal Lorentz factor (left panel) and acceleration efficiency (right panel). 


radii, before reaching the coasting regime, so that the coasting Lorentz factor (i.e., the “initial” Lorentz factor for 
the afterglow) is at least the “terminal” Lorentz factor. Here we only focus on 1-D models with Cartesian geometry. 
Implications for GRB models will be discussed in more detail when this Riemann problem is solved in conical jet 
geometry. 
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